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ABSTRACT 



We combine all available information to constrain the nature of OGLE-2005- 
BLG-071Lb, the second planet discovered by microlensing and the first in a high- 
magnification event. These include photometric and astrometric measurements 
from Hubble Space Telescope, as well as constraints from higher order effects ex- 
tracted from the ground-based light curve, such as microlens parallax, planetary 
orbital motion and finite-source effects. Our primary analysis leads to the conclu- 
sion that the host of Jovian planet OGLE-2005-BLG-071Lb is an M dwarf in the 
foreground disk with mass M = 0.46 it 0.04 Mq, distance Di = 3.2±0.4kpc, and 
thick-disk kinematics fLSR ~ 103kms~^. From the best-fit model, the planet has 
mass Mp = 3.8±0.4 Mjupiter? lies at a projected separation r± = 3.6±0.2 AU from 
its host and so has an equilibrium temperature of T ~ 55 K, i.e., similar to Nep- 
tune. A degenerate model less favored by Ax^ = 2.1 (or 2.2, depending on the 
sign of the impact parameter) gives similar planetary mass Mp = 3.4±0.4 Mjupiter 
with a smaller projected separation, r±_ = 2.1 ± 0.1 AU, and higher equilibrium 
temperature T ~ 71 K. These results from the primary analysis suggest that 
OGLE-2005-BLG-071Lb is likely to be the most massive planet yet discovered 
that is hosted by an M dwarf. However, the formation of such high-mass plane- 
tary companions in the outer regions of M-dwarf planetary systems is predicted 
to be unlikely within the core-accretion scenario. There are a number of caveats 
to this primary analysis, which assumes (based on real but limited evidence) that 
the unlensed light coincident with the source is actually due to the lens, that is, 
the planetary host. However, these caveats could mostly be resolved by a single 
astrometric measurement a few years after the event. 

Subject headings: gravitational lensing - planetary systems - Galaxy: bulge 



1. Introduction 



Microlensing provides a powerful met hod to detect extrasolar planets. Although only six 



microlens planets have been found to date 


Bond et al. 


2004 


Udalski et al. 


2005 


2006; 


Gould et al. 


2006; 


Gaudi et al. 


2008; 


Done; et al. 


I2OO9 


) , these include two 



eries. First, two of the planets are "cold Neptunes", a high discovery rate in this previously 



(Gould et al. 


2006; 


Kubas et al. 


2008) 



analog via a very high-magnification event with substantial sensitivity to multiple planets 
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indic ates that solar system analogs may be prevalent among planetary systems (iGaudi et al. 
20081 ) . Recent improvements in search techniques and fu ture major up grades should increase 
the discovery rate of microlensing planets substantially (lGaudill2008l ). 



Routine analysis of planetary microlensing light curves yields the planet/star mass ra- 
tio q and the planet-star projected separation d (in units of the angular Einstein radius). 
However, because the lens-star mass M cannot be simply extracted from the light curve, the 
planet mass Mp = qM remains, in general, equally uncertain. 

The problem of constraining the lens mass M is an old one. When microlensing ex- 
periments were initiated in the early 1990s, it was generally assumed that individual mass 
measurements would be impo ssible and tha t only statistical estimates of the lens mass scale 
could be recovered. However, iGouldl (Il992l ) pointed out that the mass and lens-source rel- 
ative parallax, tTi-ci = tt; — tt^, are simply related to two observable parameters, the angular 
Einstein radius, 9e, and the Einstein radius projected onto the plane of the observer, te. 



M 



Here, t te = A U/fR is the "microlens parallax" and k = 4G/(c^AU) 



(1) 

8.1mas/M0. See 



Gouldl (l2000bl ) for an illustrated derivation of these relations. 



In principle, 6e can be measured by comparing some structure in the light curve to a 
"standard angular ruler" on the sky. The best ex ample is light-curve distortions due to the 
finite angular radius of the source 9^ (lGouldlll994f). which usually can be estimated very well 
from its color and apparent magnitude (lYoo et al.l l2004j ). While such finite-source effects 
are rare for microlensing events considered as a whole, they are quite common for planetary 
events. The reason is simply that the planetary distortions of the light curve are typically 
of similar or smaller scale than 6^. In fact, all six planetary events discovered to date show 
such effects. Combining 6e with the (routinely measurable) Einstein radius crossing time 
yields the relative proper motion /i in the geocentric frame, 

Oe 

f^geo = (2) 
tE 

From equation ([1]), measurement of 9e by itself fixes the product Mn^ei = 6'|/k. Using priors 
on the distribution of lens-source relative parallaxes, one can then make a statistical estimate 
of the lens mass M and so the planet mass Mp. 

To do better, one must develop an additional constraint. This could be measurement 
of the microlens parallax tte, but this is typically possible only for long events. Another 
possibility is direct detection of the lens, either under the "glare" of the source during 
and immediately after the event, or displaced from the source well after the event is over. 
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Bennett et al.l (120061 ) used the latter technique to constrain the mass of the first microlensing 
planet, OGLE-2003-BLG-235/MOA-2003-BLG-53Lb. They obtained Hubble Space Tele- 
scope (HST) Advanced Camera for Surveys (ACS) images in B, V, and / at an epoch 
At = 1.78 years after the event. They found astrometric offsets of the (still overlapping) 
lens and source light among these images of up to 0.7 mas. Knowing the lens-source angular 
separation AO = /lAt from the already determined values of 9e and ts, they were able to 
use these centroid offsets to fix the color and magnitude of the lens and so (assuming that 
it was a main- sequence star) its mass. 

While the planet mass is usually considered to be the most important parameter that 
is not routinely derivable from the light curve, the same degeneracy impacts two other 
quantities as well, the distance and the transverse velocity of the lens. Knowledge of these 
quantities could help constrain the nature of the lens, that is, whether it belongs to bulge, 
the foreground disk, or possibly the thick disk or even the stellar halo. Since microlensing 
is the only method currently capable of detecting planets in populations well beyond the 
solar neighborhood, extracting such information would be quite useful. Because the mass, 
distance, and transverse velocity are all affected by a common degeneracy, constraints on 
one quantity are simultaneously constraints on the others. As mentioned above, simultane- 
ous measurements of and tte directly yield the mass. However, clearly from equation ([1]) 
they also yield the distance, and hence (from eq. [2]) also the transverse velocity. Here, we 
assemble all available data to constrain the mass, distance and transverse velocity of the sec- 
ond microlensing pl anet, OGLE-2005-BLG-071Lb, whose discovery we previously reported 



(lUdalski et al.ll2005l . hereafter Paper I). 



2. Overview of Data and Types of Constraints 

The light curve consists of 1398 data points from 9 ground-based observatories (see 
Fig. [I]), plus two epochs of HST ACS data in the F814W (/) and F555W (V) filters. The 
primary ground-based addition relative to Paper I is late-time data from OGLE, which 
continued to monitor the event down to baseline until HJD = 2453790.9. 

These data potentially provide constraints on several parameters in addition to those 
reported in Paper I. First, the light curve shows a clear asymmetry between the rising 
and falling parts of the light curve, which is a natural result of microlens parallax due to 
the Earth's accelerated motion around the Sun (see the best-fit model without parallax 
effects plotted in dotted line in Fig. [1]). However, it is important to keep in mind that such 
distortions are equall y well produced by "xall arap" due to accelerated motion of the source 



around a companion. iPoindexter et al.l (120051 ) showed that it can be difficult to distinguish 
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between the two when, as in the present case, the effect is detected at Ax^ ^100. 

Second, the two pronounced pealcs of the hght curve, which are due to "cusp approaches" 
(see the bottom inset of Fig. [T]), are relatively sharp and have good coverage. These peaks 
would tend to be "rounded out" by finite-source effects, so in principle it may be possible to 
measure p (i.e., 9^, in the units of 9e) from these distortions. 

Third, the orbital motion of the planet can give rise to two effects: rotation of the caustic 
about the center of the mass and distortion of the caustic due to expansion/contraction of 
the planet-star axis. The first changes the orientation of the caustic structure as the event 
evolves while the second changes its shape. These effects are expected to be quite subtle 
because the orbital period is expected to be of an order of 10 years while the source probes 
the caustic structure for only about 4 days. Nevertheless, they can be very important for 
the interpretation of the event. 

Finally, the HST data cover two epochs, one at 23 May 2005 (indicated by the arrow 
in Fig. [1]) when the magnification was about A = 2 and the other at 21 Feb 2006 when the 
event was very nearly at baseline, A ~ 1. These data could potentially yield four types of 
information. First, they can effectively determine whether the blended light is "associated" 
with the event or not. The blended light is composed of sources in the same photometric 
aperture as the magnified source, but that do not become magnified during the event. If 
this light is due to the lens, a companion to the lens, or a companion to the source, it should 
fall well within the ACS point spread function (PSF) of the source. On the other hand, if it 
is due to a random interloper along the line of sight, then it should be separately resolved 
by the ACS or at least give rise to a distorted PSF. Second, the HST data can greatly 
improve the estimate of the color of the blended light. The original model determined the 
source fluxes in both OGLE V and / very well, and of course the baseline fluxes are also 
quite well determined. So it would seem that the blended fluxes, which are the differences 
between these two, would also be well determined. This proves to be the case in the I band. 
However, while the source flux is derivable solely from flux differences over the light curve 



(and so is well determined from OGLE difference image analysis - DIA - IWoznialal2000l ) . the 
baseline flux depends critically on the zero point of PSF-fitting photometry, whose accuracy 
is fundamentally limited in very crowded bulge fields. The small zero-point errors turn out 
to have no practical impact for the relatively bright / background light, but are important 
for the V band. Third, one might hop e to measure a centroid shift between the two colors 



in the manner of iBennett et al.l (120061 ). Last, one can derive the source proper motion fi^ 
from HST data (at least relative to the mean motion of bulge stars). This is important, 
because the event itself yields the source-lens relative proper motion, ^tgeo- Hence, precise 
determination of requires knowledge of two proper-motion differences, first the heliocentric 
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proper motion 

f^hci = f^i- Ms, (3) 
and second, the offset between the hehocentric and geocentric proper motions 

Mhel ~ ^'■geo — j^jj ■ 

Here, is the velocity of the Earth relative to the Sun at the time of peak magnifica- 
tion to- Note that, if the lens-source relative parallax vTrei is known, even approximately, 
then the latter difference can be determined quite well, since its total magnitude is just 
0.6 masyr^^(7ri.ci/0.17mas). 



3. Constraining the Physical Properties of the Lens and its Planetary 

Companion 

In principle, all the effects summarized in § |2] could interact with each other and with 
the parameters previously determined, leading potentially to a very complex analysis. In 
fact, we will show that most effects can be treated as isolated from one another, which 
greatly facilitates the exposition. In the following sections, we will discuss the higher order 
microlens effects in the order of their impact on the ground-based light curve, starting 
with the strongest, that is, parallax effects (§ 13.11) . followed by planetary orbital motion 
(§ 13. 2p and finally the weakest, finite-source effects (§ 13.3.41) . To study these effects, we 
i mplement Markov chain Monte Carlo (MCMC) with an adaptive step-size Gaussian sampler 



( iDoran fc Mueller! |2003| ) to perform the model fitting and obtain the uncertainties of the 
parameters. The HST astrometry is consistent with no {V — I) color-dependent centroid 
shift in the first epoch, while such a shift is seen in the second epoch observations (§ 13. 4p . 
In addition, the PSF of the source shows no sign of broadening due to the blend, suggesting 
that the blend is associated with the event (§ 13. 5p . Therefore, the i/ST observations provide 
good evidence that the blend is likely due to the lens. In § 13.41 it is shown, under such an 
assumption, how the astrometry can be used in conjunction with finite-source and microlens 
parallax measurements to constrain the angular Einstein radius and proper motion (§ 13. 4p . In 
§ I3.5[ we discuss using i/ST photometric constraints in the form of penalties to the MCMC 
runs to extract the color and brightness of the blend. The results of these runs, which include 
all higher order effects of the ground-based light curve and the i/S" T photometric constraints, 
are summarized as "MCMC A" in Table [TJ Subsequently in § 13. 6^ by making the assumption 
that the blended light seen by HST is due to the lens, we combine all constraints discussed 
above to obtain physical parameters of the lens star and its planet. The corresponding best- 
fit model parameters are reported as "MCMC B" in Table [TJ The results for the physical 
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parameters from these runs are given in Table [2l Finally we discuss some caveats in the 
analysis in § 13.71 and § 13. 8[ 



3.1. Microlens Parallax Effects 



A point-source static binary-lens model has 6 "geometric- mo del" parameters: three 
"single-lens" parameters {to, Uq, t^), where we define the time of "peak" magnification 
(actually lens-source closest approach) to and the impact parameter uq with respect to the 
center of mass of the planet-star systems; and three "binary-lens" parameters {q, d, a), where 
a is the angle between the star-planet axis and the trajectory of the source relative to the 
lens. In addition, flux parameters are included to account for light coming from the source 
star (Fs) and the blend (Fb) for each dataset. In this paper, we extend the fitting by including 
microlens parallax, orbital motion and finite-source effects. Paper I reported that, within 
the context of the point-source static binary-lens models, the best-fit wide-binary {d > 1) 
solution is preferred by Ax^ = 22 over the close-binary {d < 1) solution. Remarkably, when 
we take account of parallax, finite-source and orbital effects, this advantage is no longer as 
significant. We discuss the wide/close degeneracy with more detail in § I3.6.2[ 

The microlens parallax e ffects are parara etrize d by tt f r and 7rE,N, following the geo- 
centric parallax formalism by lAn et al.l (|2002| ) and iGouldl (l2004j ) . To p roperly model the 
parallax effects, we characterize the "constant acceleration degeneracy" (jSmith et al.ll2003l ) 
by probing models with uq —uq and a —a. We find that all other parameters re- 
main essentially unchanged under this form of degeneracy. In the following sections, if not 
otherwise specified, parameters from models with positive Uq are adopted. 

As shown in Figure [21 microlens parallax is firmly detected in this event at > 8cr 
level. Not surprisingly, the error ellipse of tte is elongated toward vrE,±, i.e., the direction 
perp endicular to the position of the Sun at the peak of event, projected onto the plane of the 
sky (iGould et al.lll994j : iPoindexter et al.ll2005l ). As a result, tte^e is much better determined 
than tte n ? 



tte^e = -0.26 ± 0.05, 7rE,N = -0.30toit (wide), (5) 
tte^e = -0.27 ±0.05, tte.n = -0.36lgJ^ (close). (6) 

Xallarap (light-curve distortion from refiex motion of the source due to a binary com- 
panion) could provide an alternate explanation of the detected parallax signals. In § 13. 8[ we 
find that the best-fit xallarap parameters are consistent with those derived from the Earth's 
orbit, a result that favors the parallax interpretation. 
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3.2. Fitting Planetary Orbital Motion 

To model orbital motion, we adopt the simplest possible model, with uniform expansion 
rate b in binary separation b and uniform binary rotation rate uj. Because orbital effects are 
operative only for about 4 days, while the orbital period is of order 10 years, this is certainly 
adequate. Interestingly, the orbital motion is more strongly detected for the close solutions 
(at ^ 5.5(7 level) than the wide solutions (at ~ 3cr level), and as a result, it significantly 
lessens the previous preference of the wide solution that was found before orbital motion was 
taken into account. Further discussions on planetary orbital motion are given in § I3.6.2[ 



3.3. Finite-source effects and Other Constraints on 6e 

3.3.1. Color- Magnitude Diagram 



We follow the standard procedure to derive dereddened source color and magnitude from 
the color-magnitude diagram (CMD) of the observed field. Figure [3] shows the calibrated 
OGLE CMD (black), with the basehne source being displayed as a green point. The V — I 
color of the source can be determined in a model-independent way from linear regression of 
the /-band and V^-band observations. The /-band magnitude of the source is also precisely 
determined from the microlens model, and it is hardly affected by any higher order effects. 
The center of red clump (red) is at {V — /, /)ciump = (1-89, 15.67). The Galactic coordinates 
of the source are at (/, b) = (355.58, —3.79). Because the Galactic bulge is a bar-like structure 
that is inclined relative to the plane of the sky, the red clump d ensity at this sky position 
peaks behind the Galactic center by 0.15 mag jNishivamal bood ). Hence, we derive {V — 
-^5 -^)o,ciump = (1-00, 14.47), by adopting a Galactic distance Rq = 8kpc. We thereby obtain 
the selective and total extinction toward the source [E(y — I),Ai] = (0.89, 1.20) and thus 
Rvi = Ay / E{y — I) = 2.35. The dereddened color and magnitude of the source is {{V — 
I),I)s,o = (0.45, 18.31). From its dereddened color {V — /)o = 0.45, as well as its absolute 
magnitude (assuming it is in the bulge) Mj ~ 3.65, we conclude that the source is a main- 
sequence turnoff star. F ollowing the method o f lYoo et al.l (120041 ). we transform {V — /)o = 
0.45 to (y — K)q = 0.93 (IBessell &: Brettlll988l ). and based on the empi rical relation betwee n 
the color and surface brightness for subgiant and main-sequence stars (IKervella et al.ll2004j ). 
we obtain the angular size of the source 



9, = 0.52 X ioO-2(i9-5i-^=) ± 0.05 /xas. 



(7) 



where Ig is the apparent magnitude of the source in the / band. Other features on the CMD 
shown in Figure [3] are further discussed in § 13.51 
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3.3.2. Photometric Systematics of the Auckland Data set 

The Auckland data set's excellent coverage over the two peaks makes it particularly 
useful for probing the finite-source effects. Unlike the more drastic "caustics crossings" that 
occur in some events, the finite-source effects during "cusp approaches" are relatively subtle. 
Hence one must ensure that the photometry is not affected by systematics at the few percent 
level when determining p = 6^/6e- The Auckland photometry potentially suffers from two 
major systematic effects. 

First, the photometry of constant stars reduced by /iFUN's DoPHOT pipeline are found 
to show sudden "jumps" of up to ~ 10% when the field crossed the meridian each night. 
The signs and amplitudes of the "jumps" depend on the stars' positions on the CCD. The 
Auckland telescope was on a German equatorial mount, and hence the camera underwent a 
meridian flip. Due to scattered light, the flat-fielded images were not uniform in illumination 
for point sources, an effect that c an be corrected by making "superflats" with photometry 



of constant stars (iManfroidI Il995l ). We have constructed such "superflats" for each night 
of Auckland observations using 71 bright isolated comparison stars across the frame. The 
DoPHOT instrumental magnitude rriij for star i on frame j is modeled by the following 
equation: 

= mo,i - fixij, yij) - Zj - fwhrriij x Si (8) 

where mo,j is the corrected magnitude for star i, f{x,y) is a. biquartic illumination correction 
as a function of the (x, y) position on the CCD frame with 14 parameters, Zj is a zero- 
point parameter associated with each frame (but with Zi set to be zero), and Si is a linear 
correlation coefficient for the seeing fwhrrii j. A least-squares fit that recursively rejects 4-0" 
outliers is performed to minimize x^- The best-fit f{x,y) is dominated by the linear terms 
and has small quadratic terms, while its cubic and quartic terms are negligible. The resulting 
reduced is close to unity, and the "jumps" for all stars are effectively eliminated. We apply 
the biquartic corrections to the images and then reduce the corrected images using the DIA 
pipeline. The resulting DIA photometry of the microlens target is essentially identical (at 
the ~ 1% level) to that from the DIA reductions of the original Auckland images. 

As we now show, this is because DIA photometry automatically removes any artifacts 
produced by the first- and second-order illumination distortions if the sources are basically 
uniformly distributed across the frame. For the first-order effect, a meridian flip about the 
target (which is very close to center of the frame) will induce a change in the flux from the 
source, but it will also induce a change in the mean flux from all other stars in the frame, 
which for a linear correction will be the same as the change in position of the "center of 
light" of the frame light. If the frame sources are uniformly distributed over the frame, the 
"center of light" will be the center of the frame, which is the same position as the source. 
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therefore introducing no effects. The second-order transformation is even under a rotation 
of 180 degrees, whereas a meridian flip is odd under this transformation. Hence the flip has 
no effects at second order. 

Second, the Auckland observations were unfiltered. The amount of atmospheric extinc- 
tion differs for stars with different colors. As shown in Figure [3], the source is much bluer 
than most of the bright stars in the field, which dominate the reference image. So the amount 
of extinction for the source is different from the average extinction over the whole frame. 
This difference varies as the airmass changes over time during the observations. Coinciden- 
tally, the times of the two peaks were both near maximum airmass when the "differential 
extinction" effect is expected to be the most severe. To investigate this effect, we match the 
isolated stars in the Auckland frame with CTIO / and V photometry. We identify 33 bright, 
reasonably isolated stars with \{V — I) — {V — I)s\ < 0.25. We obtain a "light curve" for 
each of these stars, using exactly the same DIA procedure as for the source. We measure the 
mean magnitude of each of the 33 light curves and subtract this value from each of the 508 
points on each light curve, thereby obtaining residuals that are presumably due primarily to 
airmass variation. For each of the 508 epochs, we then take the mean of all of these residuals. 
We recursively remove outliers until all the remaining points are within 3 a of the mean, as 
defined by the scatter of the remaining points. Typically, 1 or 2 of the 33 points are removed 
as outliers. The deviations are well fitted by a straight line, 

dMag , , 

— = 0.0347 ±0.0016 (9) 
clZi 

where Z is airmass. The sense of the effect is that stars with the color of the microlensed 
source are systematically fainter at high airmass, as expected. (We also tried fitting the data 
to a parabola rather than a line, but the additional [quadratic] parameter was detected at 
substantially below 1 a.) Finally, we apply these "differential extinction" corrections to the 
"superflat"-adjusted DIA photometry to remove both photometric systematics. 

In general, the finite-source effects depend on the limb-darkening profile of the source 
star in the observed passbands. We find below that in this case, the impact proves to be 
extremely weak. Nevertheless, using the matched Auckland and CTIO stars, we study the 
difference between Auckland magnitudes and /-band magnitude as a function of the V — I 
color. We find the Auckland clear filter is close to the R band. 



3.3.3. Blending in Palomar and MDM Data 



Palomar data cover only about 80 minutes, but these include the cresting of the second 
peak, from which we derive essentially all information about finite-source effects. The Palo- 



- 13 - 



mar data are sensitive to these effects through their curvature. The curvature derived from 
the raw data can be arbitrarily augmented in the fit (and therefore the finite-source effects 
arbitrarily suppressed) by increasing the blending. In general, the blending at any obser- 
vatory is constrained by observations at substantially different magnifications, typically on 
different nights. However, no such constraints are available for Palomar, since observations 
were carried out on only one night. 

We therefore set the Palomar blending = 0.2 fg, that is, similar to the OGLE blend- 
ing. That is, we assume that the observed flux variation of 9%, over the Palomar night, 
actually reflects a magnification variation of 9%/[l — fb/Afg] = 9% + 0.026%, where A ~ 70 
is the approximate magnification on that night. If our estimate of the blending were in error 
by of order unity (i.e. either /(, = or /f, = 0.4 /g), then the implied error in the magni- 
fication difference would be 0.026%, which is more than an order of magnitude below the 
measurement errors. Hence, the assumption of fixed blending does not introduce "spurious 
information" into the fit even at the 1 a level. MDM data cover the second peak for only 
~ 18 minutes. For consistency, we treat its blending in the same way as Palomar, although 
the practical impact of this data set is an order of magnitude smaller. 



3.3.4- Modeling the Finite-Source Effects 



After careful tests that are described immediately below, we determined that all finite- 
source calc ulation s can b e carried out to an accuracy of 10~^ using the hexadecapole approx- 
imation of iGouldl (120081 ) (see also lPejcha fc Heyrovskyll2008l ). This sped up calculations by 
several orders of magnitude. We began by conducting MCMC sirnulatio ns using the "loop 
linking" finite-source code described in Appendix A of Dong et al. ( 20061 ). From these sim- 
ulations, we found the 4.5 cr upper bound on the finite-source parameter p[A.ba) = 0.001. 
We then examined the differences between loop-linking (set at ultra-high precision) and hex- 
adecapole for_Jight_curves at this extreme limit and found a maximum difference of 10~^. 
Based on I Clare tl (120001 ). we adopt linear limb- darkening coefficients Tj = 0.35 for the I-band 
observations and P/j = 0.43 for the observations performed in the R-band and the clear 
filters, where the local surface brightness is given by S{9) oc 1 — r[l — 1.5(1 — O'^/OlY^'^]. Ten 
additional MCMC runs are performed with P/ and Pi? that differ from the above values by 
0.1 or 0.2. They result in essentially the same probability distributions of p. Therefore, the 
choice of limb-darkening parameters has no effect on the results. The source size is found 



to be p = 3.9^*^2!? the wide solution and p 



3.lt.2l for the close solution. 



Solutions with 



p > 0.0009 are ruled out at more than 3a. The angular Einstein radius is given by 6e = 0^,/p. 
Hence, the lack of pronounced finite-source effects yields a 3a lower limit: 6e > 0.6 mas. The 
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lens-source relative proper motion in the geocentric frame is simply Hgeo = ^e/^e- The pos- 
terior probability distributions of fi^^^ derived from these MCMC simulations are compared 
with those derived from astrometry in § 13.41 



3.4. HST Astrometry 

i/5r observations were taken at two epochs (HJD = 2453513.6 and HJD = 2453788.2) 
with the ACS High Resolution Camera (HRC). For each epoch, 4 dithered images were 
acquired in each of F814W and F555W with individual exposure times of 225s and 315s 
respectively. The position of the microlens on the HST frame is in excellent agreement 
with its centroid on the OGLE difference image (within ~ 0".01). The closest star to the 
source is about 0".6 away. This implies that the OGLE photometry of the target star does 
not contain additional blended light that would be identifiable fro m the HST images. Data 



analysis was carried out using the softwar e program img2 x yni_HR C ( lAnderson fc King 12004] ) 



in a manner similar to that described in Bennett et al. (120061). Stars are fitted with an 



empirical "library" PSF that was derived from well-populated globular cluster fields. These 
positions are then corrected with precise distortion-correction models (accurate to ~ 0.01 
pixel). We adopted the first F555W frame of the first epoch as the reference frame, and 
used the measured positions of stars in this frame and the frame of each exposure to define 
a linear transformation between the exposure frame and the reference frame. This allowed 
us to transform the position of the target star in each exposure into the reference frame, so 
that we see how the target star had moved relative to the other stars. The centroid positions 
of the target star in each filter and epoch are shown in Figure |H For convenience, in this 
figure, the positions are displayed relative to the average of the centroid positions. The error 
bars are derived from the internal scatter of the four dithered images. The probability is 
P = 38% of measuring the observed separation (or larger) between F814W and F555W under 
the assumption that the true offset is zero. The fact that the blended light is aligned with the 
source argues that it is associated with the event (either it is the lens itself or a companion 
to the lens or the source). We give a more quantitative statement of this constraint in § 13.51 
For the present we simply note that the P = 38% probability is compatible with the picture 
that the blend is due to the lens since the first epoch was only about half of the Einstein- 
radius crossing time after to, implying that the lens-source separation induces only a very 
small centroid offset, well below the i/S" T detection limit. For the second epoch, the centroid 
offset is, 



^i"F8i4W-F555W,East — —0.52 ± 0.20 mas, ArF8i4W-F555W,North — 0.22 ± 0.20 mas. (10) 
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We also calculate the error in the centroid offsets from the scatter in such offsets among all 
comparison stars with F555W magnitudes within 0.5 mag of the target and find that it is 
consistent with the internally-based error quoted above. 

At the peak of the event, the angular separation between the lens and the source was 
neghgible, since mq <^ 1. We therefore fix the angular positions of the lens and source at a 
common Bq. From the CMD (Fig. [3]), most of the stars in the HST field are from the bulge. 
So we set a reference frame that is fixed with respect to the bulge field at distance Dg. The 
source and lens positions at time t are then, 



Ol{t) = Oo + ^ll{t - to) + 7T,,i[s{t) ~ s{to)], 



(11) 



where s(t) is the Earth-to-Sun vector defined by Gould (2004). Then by applying equa- 
tions ([3]) and (jl]), the angular separation between the lens and source is. 



fre^t) = Olit) - Osit) = fl„,,{t - to) + 7rrelAs(t) 



(12) 



where As(t) is given by eq. (5) in iGouldl (120041 ). 



The centroid of the source images 0'^ is displaced from the source position by (jWalker 



19951), 



Therefore, one can obtain the centroid position of the lens and the source at time t, 

Ooit) = [i-m][ds{t) + AOsit)] + moiit) 

= 6o + Mt - to) + OUtmt) + J (14) 

where fi{t) is the fraction of the total flux due to the lens. 

The centroid offset between the two passbands, F814W and F555W, is related to the 
properties of the system by, 

A0c(^)f814W-F555W = [fit) fit) «,F555W] X [1 - TT, u\ la 12 , o]^rel(t). (15) 

Frollrj/t7Ej + 2 

The difference of the blend's fractional flux between F814W and F555W is obtained from 
"MCMC A" described in § 13. 5[ Consequently, under the assumption that the blend is the 
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lens, we can use the measurement of the second-epoch HST centroid offset to estimate the 
relative proper motion from equation f|T5l) for a given tTi-ci. For purposes of illustration, we 
temporarily adopt vTrci = 0.2 when calculating the probability distribution of //gg^ (black 
contours in the upper panels of Fig. [S]). The centroid shift generally favors faster relative 
proper motion than that derived from the source size measurement (green contours in Fig. [5]), 
but the difference is only at the ~ la level. We then get a joint probability distribution of 
/Xggo from both finite-source effects and astrometry, which is shown as the red contours in 
the upper panels of Figure Ei 

We then derive the distribution of the ^^^^ position angle (p^ (North through East), 
which is shown by the red histograms in the lower panels of Figure [51 Since the direction 
of the lens-source relative proper motion ^^^^ is the same as that of the microlens parallax 
tte in the geocentric frame, we have an independent check on the 0^ from our parallax 
measurements, whose distribution is plotted as blue histograms in Figure [51 Both constraints 
favor the lens-source proper motion to be generally West, but they disagree in the North- 
South component for which both constraints are weaker. The disagreements between two 
histograms is at about 2.5cr level. 



3.5. "Seeing" the Blend with HST 

If the blend were not the lens (or otherwise associated with the event), the PSF of the 
source would likely be broadened by the blended star. We examine the i/S'rF814W images 
of the target and 45 nearby stars with simi lar brightness for e ach a vailable exposure. We 



fit them with the library PSF produced by [Anderson fc King (120041 ). In order to account 



for breathing-related changes of focus, we fit each of these 45 nearby stars with the library 
PSF, and construct a residual PSF that can be added to the library PSF to produce a PSF 
that is tailor-made for each exposure. For both epochs, the source-blend combination shows 
no detectable broadening relative to the PSFs of other isolated stars in the field. From the 
ground-based light curve, it is already known that ~ 16% of this light comes from the blend. 
We add simulated stars with the same flux as the blended light from to 2.0 pixels away from 
the center of the source. We find that the blend would have produced detectable broadening 
of the PSF if it were more than 15 mas apart from the source at the second epoch. Hence, 
the source-blend separation must then be less than about 15 mas. From the HST image 
itself, the density of ambient stars at similar magnitudes is < 1 arcsec"^. The probability of 
a chance interloper is therefore < 0.07%, implying that the blended light is almost certainly 
associated with the event, i.e., either the lens itself, a companion to the lens, or a companion 
to the source. Both of the latter options are further constrained in § 13. 71 where, in particular. 
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we essentially rule out the lens-companion scenario. 

As discussed in § [2], the blended flux in I is relatively well determined from the ground- 
based OGLE data alone, but the blended V flux is poorly determined, primarily because the 
systematic uncertainty in the zero point of the baseline flux (determined from PSF fitting) 
is of the same order as the blended flux. Because the HST image is very sparse, there is 
essentially no zero-point error in the HST V^-band flux. The problem is how to divide the 
baseline V flux into source and blend fluxes, -Fbase = Fs + Ff,. 

The standard method of doing this decomposition would be to incorporate the HST V 
light curve into the overall fit, which would automatically yield the required decomposition. 
Since this "light curve" consists of two points, the "fit" can be expressed analytically 

Fs = "^^'i'fl''^ = Fih) - F,. (16) 

where we have made the approximation that the second observation is at baseline. Let us 
then estimate the resulting errors in Fg and Ff,, ignoring for the moment that there is some 
uncertainty in Ai due to uncertainties in the general model. Each of the individual flux 
measurement is determined from 4 separate subexposures, and this permits estimates of 
the errors from the respective scatters. These are ai = 0.01 and (72 = 0.03 mag. Hence, 
the fractional error in F, is (2.5/ In 10)a(F,)/F, = [al{Ai + + (t|(1 + r)2]VV(^i - 1), 
where r = Fb/Fg. Adopting, for purposes of illustration, Ai = 2 and r = 0.1, this implies 
an error ctIVs^hst) of 0.04 mag. This may not seem very large, but after the subtraction 
in equation (fT6|) . it implies an error aiVh^HST) ~ cr(ys,HST)/r ~ 0.4 mag. And taking into 
account of the uncertainties introduced by model fitting in determining the magnifications, 
the error is expected to be even larger. Hence, we undertake an alternate approach. 

Because the HST and OGLE V filters have very nearly the same wavelength center, 
Vs,HST should be nearly identical to K,ogle up to a possible zero-point offset on their re- 
spective magnitude scales. Because the OGLE data contain many more points during the 
event, some at much higher magnification than the single HST event point, Ki,oGLE is deter- 
mined extremely well (for fixed microlensing model), much better than the 0.04 mag error for 
Vs,HST- Thus, if the zero-point offset between the two systems can be determined to better 
than 0.04 mag, this method will be superior. Although the J-band blend is much better 
measured than the V^-band blend from the ground-based data, for consistency we determine 
the zero-point offset in I by the same procedure. 

Figure E] shows differences between OGLE and HST V magnitudes for matched stars 
in the HST image. The error for each star and observatory is determined from the scatter 
among measurements of that star. We consider only points with V < 19.5 because at fainter 
magnitudes the scatter grows considerably. Each star was inspected on the HST images. 
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and those that would be significantly blended on the OGLE image were eliminated. The 
remaining points are fit to an average offset by adding a "cosmic error" in quadrature to 
the errors shown. We carry out this calculation twice, once including the "outlier" (shown 
as a filled circle) and once with this object excluded. For the V band, we find offsets of 
Vhst - VoGLE = 0.17 ± 0.01 and 0.18 ± 0.01, respectively. We adopt the following the 
V^-band offset 

= Vhst - Vqgle = 0.18 ± 0.01. (17) 

A similar analysis of the / band leads to 

AI = Ihst - /ogle = 0.08 ± 0.01. (18) 

We find no obvious color terms for either the V^-band or J-band transformations. As a check, 
we perform linear regression to compare the OGLE and HST (V^ — /) colors, and we find they 
agree within 0.01 mag, which further confirms the color terms are unlikely to be significant 
in the above transformations. 

We proceed as follows to make HST-hased MCMC ("MCMC A") estimates of Haggle 
and /fe^oGLE that place the blending star on the OGLE-based CMD. Since fiux parameters 
are linear, they are often left free and fitted by linear least-squares minimization, which 
significantly accelerates the computations. However, for "MCMC A" , the source fiuxes from 
OGLE and HST are treated as independent MCMC parameters so that they can help align 
the two photometric systems as described below. Since HST blended light is not affected by 
light from ambient stars (as OGLE is), we also leave i/S* T blended fiuxes as independent. 
Therefore, in "MCMC A", we include the following independent MCMC fiux parameters, 
Fi,s,OGLE, /^y,.,OGLE, Fi^s,HST, Fv,s,HST, Fi^b,HST, and Fv,b,HST, which for convenience we 
express here as magnitudes. For each model on the chain, we add to the light-curve based 
two additional terms Axy = {Vs,hst - K.ogle - AVf /[cr{AV)]'^ and Axj = {Is,hst - 
/s.OGLE — ^/)^/[o"(^/)]^ to enforce the measured offset between the two systems. Finally, we 
evaluate the \^-band blended flux from HST and convert it to OGLE system, Vb^ocLE/HST = 
K,OGLE — Vs^HST + Vb^HST (and similarly for / band), where all three terms on the rhs are 
the individual Monte Carlo realizations of the respective parameters. 

The result is shown in Figure [3], in which the blend (magenta) is placed on the OGLE 
CMD. Also shown, in cyan points, is HST photometry (aligned to the OGLE system) of the 
stars in the ACS subfleld of the OGLE fleld. Although this fleld is much smaller, its stars 
trace the main sequence to much fainter magnitudes. The blend falls well within the bulge 
main sequence revealed by the HST stars on the CMD, so naively the blend can be interpreted 
as being in the bulge. Hence, this diagram is, in itself, most simply explained by the blend 
being a bulge lens or a binary companion of the source. However, the measurement of V — I 
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color has relatively large uncertainty, and it is also consistent with the blend being the lens 
(or a companion to it) several kpc in front of the bulge, provided the blend is somewhat 
redder than indicated by the best-fit value of its color. In the following section, we assume 
the blended light seen by HST is the foreground lens star, and the HST photometry is 
combined with other information to put constraints on the lens star under this assumption. 



3.6. Final Physical Constraints on the Lens and Planet 

3.6.1. Constraints on a Luminous Lens 

In the foregoing, we have discussed two types of constraints on the host star properties: 
the first class of constraints, consisting of independent measurements of tte, and /it, relate 
the microlens parameters to the physical parameters of the lens; the second class are HST 
and ground-based observations that determine the photometric properties of the blend. 

In this section, we first describe a new set of MCMC simulations taking all these con- 
straints into account. Similarly to what is done to include i/ST photometry in the "MCMC 
A" (see § 13.51) . we incorporate HST astrometry constraints by adding penalties to the 
fittings. For a given set of microlens parameters, we can derive the physical parameters, 
namely, M, vTrei, Mgec and so calculate p = 6^./6-e (from eq. [1]) and the F814W — F555W 
centroid offset (from eq. [E]). Then we assign the penalties based on the observed centroid 
offset from § 13. 4[ In this way, the MCMC simulations simultaneously include all microlens 
constraints on the lens properties. The posterior probability distribution of M and vTrei are 
plotted in Figure [71 The vTrei determination very strongly excludes a bulge (vTrei < 0.05) lens. 
Note that by incorporating HST astrometry, we implicitly assume that the blend is the lens. 

If the blend is indeed the lens itself, we can also estimate its mass and distance from the 
measured color and magnitude of the blend. In doing so, we use theoretical stellar isochrones 
( M. Pinsonneault 2007, priva te communication) incorporating the color-temperature relation 



by lLejeune et al.l (119971 . Il998l ). We first use an isochrone that has solar metal abundance, with 
stellar masses ranging from O.25M0 — I.OMq, and an age of 4Gyrs. The variation in stellar 
brightness due to stellar age is negligible for our purpose. Extinction is modeled as a function 
of Di by dAj/dDi = (0.4kpc~^) exp(— wZ);), where w is set to be 0.31kpc~^ so that the 
observed value A/(8.6kpc) = 1.20 (as derived from CMD discussed in § l3.3.ip is reproduced. 
Again, the distance to the source is assumed to be 8.6 kpc, implying vr^ = 0.116 mas, and 
hence that the lens distance is Di/kpc = mas/ (vTrei + vr^). In Figure [3, we show the lens mass 
M and relative parallax vTrei derived from the isochrone that correspond to the observed I- 
band magnitude / = 21.3 in black line and a series of V — I values V — I = 1.8 (best estimate). 
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2.0 (0.5 cr), 2.1 (1 a), 2.3 (1.5 cr) and 2.6 (2 a) as black points. The observed color is in modest 
disagreement < 2cr with the mass and distance of the lens at solar metallicity. We also 
show analogous trajectories for [M/H] = —0.5 (red) and [M/H] = —1.0 (green). The level of 
agreement changes only very weakly with metallicity. 

We then include the isochrone information in a new set of MCMC runs ("MCMC B"). 
To do so, the HST blended fluxes in / and V bands can no longer be treated as independent 
MCMC parameters. Instead, based on the isochrone with solar metallicity, the lens V ~ I 
color and I magnitude are predicted at the lens mass and distance determined from MCMC 
parameters. Then the HST /-band and l^-band fluxes are fixed at the predicted values in 
the fitting for each MCMC realization. 

Figure [S] illustrates the constraints on M and vTrei from the MCMC, which are essentially 
the same for both wide-binary (solid contours) and close-binary (dashed contours) solutions: 

M = 0.46 ± 0.04 M0, TTrei = 0.19 ±0.03 mas. (19) 

Assuming the source distance at 8.6 kpc, the n^ei estimates translate to the following lens 
distance measurement: 

A = 3.2 ±0.4 kpc. (20) 

Furthermore, we can derive constraints on the planet mass Mp and the projected sepa- 
ration between the planet and the lens star r±, 

Mp = 3.8 ± 0.4Mjupiter, r± = 3.6 ± 0.2AU (wide), (21) 

and 

Mp = 3.4 ± 0.4Mjupiter, = 2.1 ± O.lAU (close). (22) 

The wide solution is slightly preferred over close solution by Ax^ = 2.1. 

To examine possible uncertainties in extinction estimates, we reran our MCMC with Aj 
and Av that are 10% higher and lower than the fiducial values. These runs result in very 
similar estimates as when adopting the fiducial values. 

From equations (IT^ and (1131) . one can easily obtain the centroid shift between two 
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epochs in a given passband by ignoring A.Os{t) 0, 

Oa{t2)-0^{h) = fJ,s{t2-tl)+tl^Jfl{t2){t2-to)- fl{h){ti-to)] 

+7r,ei[/Kt2)As(t2) - /i(tl)As(ti)] (23) 

Because /Xgg^, vTrei and fi in a given passband can be extracted from the MCMC reahza- 
tions ("MCMC B"), we can use the above equation to measure the source proper motion by 
making use of the centroid shift in F814W between two epochs. The source proper motion 
with respect to the mean motion of stars in the HST field is measured to be 

t^s = {f^s,E, fis,N) = (2.0 ± 0.2, -0.51°:?) masyr-i. (24) 

We obtain similar results with F555W, but with understandably larger errorbars since the 
astrometry is more precise for the microlens in F814W. 



Combining equations ([3]) and (jl]), the lens proper motion in the heliocentric frame is 
therefore 

fJ-i = Mgco + Ms + (25) 

For each MCMC realization, vTrei is known, so we can convert the lens proper motion to the 
velocity of the lens in the heliocentric frame v^^hei and also in the frame of local standard 
of rest Vi,LSR (we ignore the rotation of the galactic bulge). The lens velocity in the LSR is 
estimated to be f;,LSR = 103 ± 15kms~^. This raises the possibility of the lens being in the 
thick disk, in which the stars are typically metal-poor. As shown in Figure [TJ the constraints 
we have cannot resolve the metallicity of the lens star. 



3.6.2. Planetary Orbital Motion 



Wide/Close De generacy Binary- lens li ght curves in general exhibit a well-known "close- 
wide" symm etry ( DominikI 1 1 999l : lAnI l2005l ). Even for some well-covered caustics-crossing 
events (e.g., Albrow et al.l 1999 ). there are quite degenerate sets of solutions between wide 



and close binaries. In Paper I, we found that the best-fit point-source wide-binary solution 



^ The angular separations between the source and the lens are ~ O.470E and ~ 4.46'e for the two HST 
epochs, respectively. Thus the angular position offsets between the centroids of the source images and the 
source are both ^ 0.216'e and the directions of the offset relative to the source are almost the same due to 
the small impact parameter uq. The difference between lens flux fractions of the two epochs are about 7% 
in / band, so the offsets can be confidently ignored in deriving the source proper motion using the relative 
astrometry in F814W at two different epochs. 
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was preferred over close-binary solutions by Ax^ ~ 22. But this did not necessarily mean 
that the wide-close binary degeneracy was broken, since the two classes of binaries may be 
influenced differently by higher order effects. We find that the difference between best-fit 
wide and close solutions is within 1 from "MCMC A" and 2.1 (positive Uq) or 2.2 (negative 
Mo) from "MCMC B". 

However, orbital motion of the planet is subject to additional dynamical constraints: 
the projected velocity of the planet should be no greater than the escape velocity of the 
system: v±_ < v^sc, where, 

I ATT 

= Jd^ + [ujdf ^^E, (26) 



2GM 2GM ri^ , , 

and where r is the instantaneous 3-dimensional planet-star physical separation. Note that 
in the last step, we have used equation ([1]). 



We then calculate the probability distribution of the ratio 



(28) 



c2 [tte + (vr,/0E)]=^ ^^E 

for an ensemble of MCMC realizations for both wide and close solutions. Figure [9] shows 
probability distributions of the projected velocity in the units of critical velocity fc,±, 
where r_|_7 is the instantaneous velocity of the planet on the sky, which is further discussed in 
Appendix |X] and fc,± = fesc,±/v^- The dotted circle encloses the solutions that are allowed 
by the escape velocity criteria, and the solutions that are inside the solid line are consistent 
with circular orbital motion. We find that the best-fit close-binary solutions are physically 
allowed while the best-fit wide-binary solutions are excluded by these physical constraints 
at 1.6(7. The physically excluded best-fit wide solutions are favored by Ax^ = 2.1 (or 2.2) 
over the close solutions, so by putting physical constraints, the degenerate solutions are 
statistically not distinguishable at 1 a. 



Circular Planetary Orbits and Planetary Parameters Planetary deviations in mi- 
cro lensing light curves are intrinsically short, so in most cases, only the instantaneous pro- 
jected distance between the planet and the host star can be extracted. As shown in § 13.6.21 for 
this event, we tentatively measure the instantaneous projected velocity of the planet thanks 
to the relatively long (~ 4 days) duration of the planetary signal. One cannot solve for the 
full set of orbital parameters just from the instantaneous projected position and velocity. 
However, as we show in Appendix |Al we can tentatively derive orbital parameters by assum- 
ing that the planet follows a circular orbit around the host star. In Figure [TOl we show the 
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probability distributions of the semimajor axis, inclination, amplitude of radial velocity, and 
equilibrium temperature of the planet derived from "MCMC B" for both wide and close so- 
lutions. The equilibrium temperature is defined to be Tgq = (Lboi/ -^^boi,©) 
where Lboi is the bolometric luminosity of the host, a is the planet semimajor axis, and i^boi,©, 
i?0, and Tq are the luminosity, radius, and effective temperature of the Sun, respectively. 
This would give the Earth an equilibrium temperature of Teq = 285 K. In calculating these 
probabilities, we assign a flat (Opik's Law) prior for the semimajor axis and assume that the 
orbits are randomly oriented, that is, with a uniform prior on cosi. 

3.7. Constraints on a Non-Luminous Lens 

In § 13.51 we noted that the blended light must lie within 15 mas of the source: otherwise 
the i/ST images would appear extended. We argued that the blended light must be associated 
with the event (either the lens itself or a companion to either the source or lens), since 
the chance of such an alignment by a random field star is < 0.07%. In fact, even stronger 
constraints can be placed on the blend-source separation using the arguments of § 13. 4[ These 
are somewhat more complicated and depend on the blend-source relative parallax, so we do 
not consider the general case (which would only be of interest to further reduce the already 
very low probability of a random interloper) but restrict attention to companions of the 
source and lens. We begin with the simpler source-companion case. 

3.7.1. Blend As Source Companion 

As we reported in § 13. 4[ there were two HST measurements of the astrometric offset 
between the V and I light centroids, dating from 0.09 and 0.84 years after peak, respectively. 
In that section, we examined the implications of these measurements under the hypothesis 
that the blend is the lens. We therefore ignored the first measurement because the lens source 
separation at that epoch is much better constrained by the microlensing event itself than by 
the astrometric measurement. However, as we now examine the hypothesis that the blend 
is a companion to the source, both epochs must be considered equally. Most of the weight 
(86%) comes from the second observation, partly because the astrometric errors are slightly 
smaller, but mainly because the blend contributes about twice the fractional light, which 
itself reduces the error on the inferred separation by a factor of 2. Under this hypothesis, we 
find a best-fit source-companion separation of 5 mas, with a companion position angle (north 
through east) of 280°. The (isotropic) error is 3 mas. Approximating the companion-source 
relative motion as rectilinear, this measurement strictly applies to an epoch 0.73 years after 
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the event, but of course the intrinsic source-companion relative motion must be very small 
compared to the errors in this measurement. 

There would be nothing unusual about such a source-companion projected separation, 
roughly 40 ± 25 AU in p hysical units. Indeed, the lo cal G-star binary distribution function 



peaks close to this value (IDuquennoy fc Mayorlll99ll ). 



The derived separation is also marginally consistent with the companion generating a 
xallarap signal that mimics the parallax signal in our dominant interpretation. The semi- 
major axis of the orbit would have to be about 0.8 AU to mimic the 1-year period of the 
Earth, which corresponds to a maximum angular separation of about 100 yuas, which is 
compatible with the astrometric measurements at the 1.6 a level. 

Another potential constraint comes from comparing the color difference with the mag- 
nitude difference of the source and blend. We find that the source is about 0.5 ± 0.5 mag too 
bright to be on the same main sequence. However, first, this is only a 1 a difference, which is 
not significant. Second, both the sign and magnitude of the difference are compatible with 
the source being a slightly evolved turnoff star, which is consistent with its color. 

The only present evidence against the source-companion hypothesis is that the astro- 
metric offset between V and / HST images changes between the two epochs, and that the 
direction and amplitude of this change is consistent with other evidence of the proper mo- 
tion of the lens. Since this is only a P = 1.7% effect, it cannot be regarded as conclusive. 
However, additional i/ST observations at a later epoch could definitively confirm or rule out 
this hypothesis. 



3.7.2. Blend As A Lens Companion 

A similar, but somewhat more complicated line of reasoning essentially rules out the 
hypothesis that the blend is a companion to the lens, at least if the lens is luminous. The 
primary difference is that the event itself places very strong lower limits on how close a 
companion can be to the lens. 



A companion with separation (in units of ^e) d ^ 1 induces a IChang-Refsdall (119791 ) 
caustic, which is fully characterized by the gravitational shear 7 = q/d^. We find that the 
light-curve distortions induced by this shear would be easily noticed unless 7 < 0.0035, that 
is, 

7 = I = ^ < 0-0035, (29) 
where Qc = M^/M is the ratio of the companion mass to the lens mass and d^. = 6c/ 6e is the 
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ratio of the lens-companion separation to the Einstein radius. Equivalently, 



V1.3 



) 



1/2 



(30) 



Here, we have normahzed qc to the minimum mass ratio required for the companion to 
dominate the hght assuming that both are main-sequence stars. (We will also consider 
completely dark lenses below). 

We now show that equation (!30!) is inconsistent with the astrometric data. If a lens 
companion is assumed to generate the blend light, then essentially the same line of reasoning 
given in § 13.7.11 implies that 0.73 years after the event, this companion lies 5 mas from the 
source, at position angle 280° and with an isotropic error of 3 mas. The one wrinkle is that 
we should now take account of the relative-parallax term in equation f|T5|) . whereas this was 
identically zero (and so was ignored) for the source- companion case. However, this term is 
only about l.SvTrei and hence is quite small compared to the measurement errors for typical 
TTrci ^ 0.2mas. We will therefore ignore this term in the interest of simplicity, except when 
we explicitly consider the case of large vTrei further below. 

Of course, the lens itself moves during this interval. From the parallax measurement 
alone (i.e. without attributing the V/ 1 astrometric displacement to lens motion), it is known 
that the lens is moving in the same general direction, i.e., with position angle roughly 210°. 
In assessing the amplitude of this motion we consider only the constraints from finite-source 
effects (and ignore the astrometric displacement). These constraints yield a hard lower limit 
on (from lack of pronounced finite-source effects) of > 0.6 mas, which corresponds to 
a proper motion yU = 3.1 mas yr^^. At this extreme value (and allowing for 2 a uncertainty 
in the direction of lens motion as well in the measurement of the companion position), the 
maximum lens-companion separation is 11.4 mas (i.e., 19 ^e), which is just ruled out by 
equation (!30|) . At larger ^e, the lens-companion scenario is excluded more robustly. For 
example, in the limit of large 6'e, we have 6'c = /i x 0.73 yr = 6'e(0.73 yr/^E) = 3.96'e, which 
is clearly ruled out by equation (1501) . 

Then we note that any scenario involving values of TTrei that are large enough that they 
cannot be ignored in this analysis (vTrei > 0.5 mas), must also have very large 6'e = t^vcx/t^'e ^ 
1 mas, a regime in which the lens-companion is easily excluded. 

The one major loophole to this argument is that the lens may be a stellar remnant 
(white dwarf, neutron star, or black hole), in which case it could be more massive than the 
companion despite the latter's greater luminosity. 
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3.8. Xallarap Effects and Binary Source 

Binary source motion can give rise to distortions of tlie liglit curve, called "xallarap" 
effects. One can always find a set of x allarap parameters to perfectly mimic parallax dis- 



tortions caused by the Earth's motion (jSmith et al.ll2003l ). However, it is a priori unlikely 



for the binary source to have such parameters, so if the parallax signal is real, one would 
expect the xallarap fits to converge to the Earth parameters. For simplicity, we assume that 
the binary source is in circular orbit. We extensively search the parameter space on a grid 
of 5 xallarap parameters, namely, the period of binary motion P, the phase A and comple- 
ment of inclination j3 of the binary orbit, which corresponds to the ecliptic longitude and 
latitude in the parallax interpretation of the light curve, as well as (^e,e, Ce,n), which are the 
counterparts of (tte.e, tte.n) of the microlens p arallax. We take advantage of the two exact 



degeneracies found by iPoindexter et al.l (120051 ) to reduce the range of the parameter search. 



One exact degeneracy takes A' = A -|- tt and Xe = ~Xe^ while all other parameters remain 
the same. The other takes (3' = — /3, Uq = —Uq and Ce n — ~^e,n (the sign of a should be 
changed accordingly as well). Therefore we restrict our search to solutions with positive Uq 
and with vr < A < 27r. In modeling xallarap, planetary orbital motion is neglected. In Fig- 
ure [HI the distribution for best-fit xallarap solutions as a function of period is displayed 
in a dotted line, and the xallarap solution with a period of 1 year has a Ax^ = 0.5 larger 
than the best fit at 0.9 year. Figure [T2] shows that, for the xallarap solutions with period 
of 1 year, the best fit has a Ax^ = 3.2 less than the best-fit parallax solution (displayed as 
a black circle point) and its orbital parameters are close to the ecliptic coordinates of event 
(A = 268°, j3 = —11°). Therefore, the overall best-fit xallarap solution has Ax^ = 3.7 smaller 
than that of the parallax solution (whose value is displayed as a filled dot in Fig. [TTl) for 
3 extra degrees of freedom, which gives a probability of 30%. The close proximity between 
the best-fit xallarap parameters and those of the Earth can be regarded as good evidence 
of the parallax interpretation. The slight preference of xallarap could simply be statistical 
flu ctuation or reflect low-l evel systematics in the light curve (commonly found in the analysis 



by IPoindexter et al.ll2005l ). 



We also devise another test on the plausibility of xallarap. In § 13.51 we argued that the 
blend is unlikely to be a random interloper unrelated to either the source or the lens. If the 
source were in a binary, then the blend would naturally be explained as the companion of the 
source star. Then from the blend's position on the CMD, its mass would be rric ~ 0.9 Mq. 
By definition, C,e is the size of the source's orbit in the units of (the Einstein radius 
projected on the source plane), 

tts arric , , 

where a is the semimajor axis of the binary orbit, and and nic are the masses of the 
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source and its companion, respectively. Then we apply Kepler's Third Law: 



2 m'^ 

(32) 



P\2 ml _ f^ErE\ 



yr7 MqItjIc + V AU / 

Once the masses of the source and companion are known, the product of and te are 
determined for a given binary orbital period P. And in the present case, fE/AU = 6eDs = 
= 4.5 X 10~^/p. By adopting = IMq, nic = 0.9Mq, for each set of P and p, 
there is a uniquely determined ^e from equation (l32i) . We then apply this constraint in the 
xallarap fitting for a series of periods. The minimum x^s for each period from the fittings 
are shown in solid line in Figure [TTl The best-fit solution has Ax^ ~ 1.0 less than the 
best-fit parallax solution for two extra degrees of freedom. Although as compared to the 
test described in the previous paragraph, the current test implies a higher probability that 
the data are explained by parallax (rather than xallarap) effects, it still does not rule out 
xallarap. 



4. Summary and Future Prospects 

Our primary interpretation of the OGLE-2005-BLG-071 data assumes that the light- 
curve distortions are due to parallax rather than xallarap and that the blended light is due 
to the lens itself rather than a companion to the source. Under these assumptions, the lens 
is fairly tightly constrained to be a foreground M dwarf, with mass M = 0.46 ± 0.04 M© 
and distance Di = 3.2 ± 0.4kpc, which has thick-disk kinematics (vlsr ~ 103kms~^). As 
we discuss below, future observations might help to constrain its metallicity. The microlens 
modeling suffers from a well-known wide-close binary degeneracy. The best-fit wide-binary 
solutions are slightly favored over the close-binary solutions, however, from dynamical con- 
straints on planetary orbital motion, the physically allowed solutions are not distinguishable 
within la. For the wide-binary model, we obtain a planet of mass Mp = 3.8 ± 0.4Mjupiter 
at projected separation rj_ = 3.6 ± 0.2 AU. The planet then has an equilibrium temperature 
of about T = 55 K, i.e. similar to Neptune. In the degenerate close-binary solutions, the 
planet is closer to the star and so hotter, and the estimates are: Mp = 3.4 ± 0.4Mjupiter, 
= 2.1 ± 0.1 AU and T ~ 71 K. 

As we have explored in considerable detail, it is possible that one or both of these 
assumptions is incorrect. However, future astrometric measurements that are made after 
the lens and source have had a chance to separate, will largely resolve both ambiguities. 
Moreover, such measurements will put much tighter constraints on the metallicity of the 
lens (assuming that it proves to be the blended light). 
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First, the astrometric measurements made 0.84 yr after the event detected motion sug- 
gests that there was still 1.7% chance that the blend did not move relative to the source. 
A later measurement that detected this motion at higher confidence would rule out the hy- 
pothesis that the blend is a companion to the source. We argued in § 13.7.11 that the blend 
could not be a companion to a main-sequence lens. Therefore, the only possibilities that 
would remain are that the lens is the blend, that the lens is a remnant (e.g., white dwarf), or 
that the blend is a random interloper (probability < 10^^). As we briefiy summarize below, 
a future astrometric measurement could strongly constrain the remnant-lens hypothesis as 
well. 

Of course, it is also possible that future astrometry will reveal that the blend does not 
moving with respect to the source, in which case the blend would be a companion to the 
source. Thus, either way, these measurements would largely resolve the nature of the blended 
hght. 

Second, by identifying the nature of blend, these measurements will largely, but not 
entirely, resolve the issue of parallax vs. xallarap. If the blend proves not to be associated 
with the source, then any xallarap-inducing companion would have to be considerably less 
luminous, and so (unless it were a neutron star) less massive than the rric = 0.9 Mq that we 
assumed in evaluating equation (132|) . Moreover, stronger constraints on te (rhs of eq. [32] ) 
would be available from the astrometric measurements. Hence, the xallarap option would 
be either excluded or very strongly constrained by this test. 

On the other hand, if the blend were confirmed to be a source companion, then essen- 
tially all higher order constraints on the nature of the lens would disappear. The parallax 
"measurement" would then very plausibly be explained by xallarap, while the "extra informa- 
tion" about 6e that is presently assumed to come from the blend proper-motion measurement 
would likewise evaporate. 

These considerations strongly argue for making a future high-precision astrometric mea- 
surement. Recall that in the il/^ST measurements reported in § I3.5[ the source and blend were 
not separately resolved: the relative motion was inferred from the offset between the V and 
/ centroids, which are displaced because the source and blend have different colors. Due 
to its well-controlled PSF, HST is capable of detecting the broadening of the PSF even if 
the separation of the lens and source is a fraction of the FWHM. Assuming that the proper 
motion is figeo ~ 4.4masyr~^, and based on our simulations in ^ 13.51 such broadening would 



be confidently detectable about 5 years after the event (see also lBennett et al.l 120071 for an- 
alytic PSF broadening estimates). Ten years after the event, the net displacement would 
be ~ 40 mas. This compares to a diffraction limited FWHM of 40 mas for H band on a 
ground-based 10m telescope and would therefore enable full resolution. The I — H color of 
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the source is extremely well determined (0.01 mag) from simultaneous / and H data taken 
during the event from the CTIO/SMARTS 1.3m in Chile. Hence, the flux allocation of the 
partially or fully resolved blend and source stars would be known. The direct detection of a 
part ially or fully resolved l ens will provide precise photometric and astrometry measurements 



see 



Kozlowski et al.l 120071 for one such example), which will enable much tighter constraints 
on the mass, distance and projected velocity of the lens. It also opens up the possibility of 
determining the metallicity of the host star by taking into account both non-photometric 
and photometric constraints. If, as indicated by the projected velocity measurem ent, it is a 



thick-disk star, then it will be one of the few such stars found to harbor a planet (iHaywood 



20081 ). 



As remarked above, a definitive detection of the blend's proper motion would still leave 
open the possibility that it was a companion to the lens, and not the lens itself. In this 
case, the lens would have to be a remnant. Without going into detail, the astrometric 
measurement would simultaneously improve the blend color measurement as well as giving 
a proper-motion estimate (albeit with large errors because the blend-source offset at the 
peak of the event would then not be known). It could then be asked whether the parallax, 
proper-motion, and photometric data could be consistently explained by any combination of 
remnant lens and main-sequence companion. This analysis would depend critically on the 
values of the measurements, so we do not explore it further here. We simply note that this 
scenario could also be strongly constrained by future astrometry. 



5. Discussion 

With the measurements presented here, and the precision with which these measure- 
ments allow us to determine the properties of the planet OGLE-2005-BLG-071Lb and its 
host, it is now possible to place this system in the context of similar planetary systems 
discovered by radial velocity (RV) surveys. Of course, the kind of information that can be 
inferred about the planetary systems discovered via RV differs somewhat from that presented 
here. For example, for planets discovered via RV, it is generally only possible to infer a lower 
limit to the planet mass, unless the planets happen to transit or produce a detectable astro- 
metric signal. Mutatis mutandis, for planets discovered via microlensing, it is generally only 
possible to measure the projected separation at the time of the event, even in the case for 



which the microlensing mass degeneracy is broken as it is here (although see iGaudi et al. 



20081 ). 



With these caveats in mind, we can compare the properties of OGLE-2005-BLG-071Lb 
and its host star with similar RV systems. It is interesting to note that the fractional 
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uncertainties in the host mass and distance of OGLE-2005-BLG-071Lb are comparable to 
those of some of the systems hsted in Table [31 

OGLE-2005-BLG-071Lb is one of only eight Jovian-mass (0.2Mjupiter < Mp < 13M.T„pi,p/ 



planets that have been detected orbiting M dwarf hosts (i.e., < 0.55 M^l (IMarcy et al. 



19981 . 12OOII : iDelfosse et al.lll998l : iButler et al.ll2006l : Ijohnson et al.ll2007bl : iBailev et al.ll2008h . 
Table [3] summarizes the planetary and host-star properties of the known M dwarf/ Jovian- 
mass planetary systems. OGLE-2005-BLG-071Lb is likely the most massive known planet 
orbiting an M dwarf. 

As suggested by the small number of systems listed in Tabled and shown quantitatively 
by several recent studies, the frequency of relatively short-period P < 2000 days, Jupiter- 
mass companion s to M dwarfs appears to be ~ 3 — 5 times lower than s uch companions 



toFGK dwarfs flButler et al.ll2006l : lEndl et al.ll2006l : I Johnson et al.ll2007bl : ICumming et al. 



20081 ) . This paucity, which has been shown to be statistically significant, is expected in the 
core-accretion model of planet formation, which generally predicts that Jovian companions 
to M dwarfs should be rare, since for lower mass stars, the dynamical time at the sites 
of planet formatio n is longer, whereas the amount of raw material available for pl anet for- 
mation is sraaller (ILaughlin et al.l |2004J : llda &: LinI l2005l : Kennedy &: KenyonI |2008| . but see 
Kornet et al.ll2006l ). Thus, these planets typically do not reach sufficient mass to accrete 
a massive gaseous envelope over the lifetime of the disk. Consequently, such models also 
predict that in the outer regions of their planetary systems, lower m ass stars should host 



a much larger p opulation of 'failed Jupiters,' cores of mass < 10 (ILaughlin et al.l 12004 



Ida fc Linll2005r). Such a population was indeed ide ntified based on two microlensing planet 
discoveries (IBeaulieu et al.ll2006l : iGould et al.ll2006l ). 



Our detection of a ~ 4Mjupiter companion to an M dwarf may therefore present a 
difficulty for the core-accretion scenario. While we do not have a constraint on the metal- 
licity of the host, the fact that it is likely a member of the thick disk suggests that its 
metallicity may be subsolar. If so, this would pose an additional difficulty for the core- 
accretion s cenario, which a lso predicts that massive planets should be r arer around raetal- 



poor stars (llda fc Linll2004[) . as has been demonstrated observationally (ISantos et al.l 12004 



Fischer &: Valentil |2005| ) . This might imply that a different mechanism is responsible for 
planet form ation in the OGL E-2005-BLG-071L system, such as the gravitational instability 
mechanism jBosslbood l2006h . 



One way to escape these potential difficulties is if the host lens is actually a stel- 
lar remnant, such as white dwarf. T he progenitors o f remnants are generally rn ore mas- 
sive sta rs, which are both p redicted (llda &: LinI l2005l : Kennedy &: KenyonI |2008| ) and ob- 
served ( Johnson et al. 2007a| jbl) to have a higher incidence of massive planets. As we dis- 
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cussed above, future astrometric measurements could constrain both the low-metaUicity and 
remnant-lens hypotheses. These measurements are therefore critical. 

Although it is difficult to draw robust conclusions from a single system, the re are now 
four published detections of Jovian-mass planetary companions with microlensing (IBond et al. 
2004 : iGaudi et al.ll2008l ). and several additional such planets have been detected that are 
currently bei ng analyzed. It is therefore reasonable to expect several detections per year 
(IGouldl 120091 ). and thus that it will soon be possible to use microlensing to constrain the 
frequency of massive planetary companions. These constraints are complementary to those 
from RV, siiice the microlensing detection method is less biased with respect to host star 
mass (lGouldll2000a] ). and furthermore probes a different region of parameter space, namely 
cool planets beyond the snow line with equilibri um temperatures similar to the giant planets 
in our solar system (see, e.g. iGould et al.l 120071 and lGouldll2009l ). 
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A. Extracting Orbital Parameters for Circular Planetary Orbit 

OGLE-2005-BLG-071 is the first planetary microlensing event for which the effects of 
planetary orbital motion in the light curve have been fully analyzed. The distortions of 
the light curve due to the orbital motion are modeled by uj and h as discussed in § 13.21 
In addition, the lens mass M and distance Di are determined, so we can directly convert 
the microlens light-curve parameters that are normalized to the Einstein radius to physical 
parameters. In this section, we show that under the assumption of a circular planetary 
orbit, the planetary orbital parameters can be deduced from the light-curve parameters. Let 
= DiO^d be the projected star-planet separation and let be the instantaneous planet 
velocity in the plane of the sky, i.e. r_|_7_|_ = rj_uj is the velocity perpendicular to this axis and 
r_L7|| = r±d/d is the velocity parallel to this axis. Let a be the semi-major axis and define 
the i, j, k directions as the instantaneous star-planet-axis on the sky plane, the direction 
into the sky, and k = i x j. Then the instantaneous velocity of the planet is 



[cos 9k + sin 9 {cos (f)i — sin (f)j)], (Al) 

where (f) is the angle between star-planet-observer (i.e., r± = a sin (p) and 9 is the angle of the 
velocity relative to the k direction on the plane that is perpendicular to the planet-star-axis. 
We thus obtain 

[GMcos9 [GM . ^ ^ 

1± = \ — 711 = \ ^^sm 9 cot (p. (A2) 
V a-^ sm0 V a-^ 

To facilitate the derivation, we define 

A = ^ = -tan^cos0, B = ^^^^ = cos^^sin^, (A3) 
7j_ GM 

which yield as an equation for sin 0: 

5 = F(sin0); ^(^) = J^^^- (A4) 



Note that F'(sin0) = when sin^ 0, = (3/2)^^ + 1 - \A\^/{9/4:)A^ + 2 . So equation M 
has two degenerate solutions when B < F(sin0^,) and has no solutions when B > F(sin0^,). 
Subsequently, one obtains. 



a = , cosz = — sm0cosy, K = \ qsmi, (,A5j 

sin V a 

where i is the inclination and K is the amplitude of radial velocity. 



The Jacobian matrix used to transform from P(r_|_, 7_|_, 7||) to P(a, 0, 9) is given below. 



-33- 



^(r±,7±,7||) 
d{a, 0, 9) 



GM 



GM 



sin 

3 cose 
2a sin </> 

— sin^cot0 



cot' 



acos(f) 

_ cos cos (, 
sin^ <l> 
sin 9 
sin^ <f> 

,2 /) j.„„2 . 



COS 9 cot 



sin'' 9 tan'' 



')■ 



Then for an arbitrary function H{a), 

5(rx,7±,7||) 5(?^±,7±,7||) 



X 



d{H{a),cos(f),9) d{a,(l),9) sin (j)H' (a)' 
which, for the special case of a fiat distribution, H{a) — In a, yields. 



'9(r^,7±,7||) 
d{ln{a), cos 4>, 9) 



GM cos2 . 



sm ( 



sin ^tan' 



(A6) 



(A7) 



(A8) 
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OGLE-2005-BLG-071 




3440 3460 3480 3500 3520 
HJD' = HJD - 2450000 

Fig. 1. — Main panel: all available ground-based data of the microlensing event OGLE-2005- 
BLG-071. HST ACS HRC observations in F814W and F555W were taken at two epochs, 
once when the source was magnified by A ~ 2 [arrow), and again at HJD = 2453788.2 (at 
baseline). Planetary models that include (solid) and excludes (dotted) microlens parallax are 
shown. Zoom at bottom: triple-peak feature that reveals the presence of the planet. Each 
of the three peaks corresponds to the source passing by a cusp of the central caustic induced 
by the planet. Upper inset: trajectory of the source relative to the lens system in the units 
of angular Einstein radius 6e- The lens star is at (0,0), and the star-planet axis is parallel 
to the X-axis. The best-fit angular size of the source star in units of is p ~ 0.0006, too 
small to be resolved in this figure. 



- 38 - 



- 



o 



0.5 - 



-1 








-0.5 



7T 



E.East 



Fig. 2. — Probability contours (Ax^ = 1,4) of microlens parallax parameters derived 
from MCMC simulations for wid e-binary (in s olid line) and close-binary (in a dashed line) 



solutions. Fig. 2 and eq. (12) in 
nE,± form a right-handed coordinate system. 



Gouldl (12004 ) imply that 7!'e,± is defined so that nE,\\ and 
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Fig. 3.— CMD for the OGLE-2005-BLG-071 field. Black dots are the stars with the OGLE 
/-band and \^-band observations. The red point and green points show the center of red 
clump and the source, respectively. The errors in their fiuxes and colors are too small to be 
visible on the graph. Cyan points are the stars in the ACS field, which are photometrically 
aligned with OGLE stars using 10 common stars. The magenta point with error bars show 
the color and magnitude of the blended light. 
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Fig. 4. — HST ACS astrometric measurements of the target star in F814W (red) and F555W 
(blue) filters in 2005 (filled dots) and 2006 (open dots). The center positions of the big circles 
show mean values of the 4 dithered observations in each filter at each epoch while radii of 
the circle represent the 1 a errors. 
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Fig. 5. — Upper two panels sliow posterior probability contours at Ax^ = 1 (solid line) and 
4 (dotted line) for relative lens-source proper motion Hg^o- The left panel is for wide-binary 
solutions and the right one is for close-binary. The green contours show the probability 
distributions constrained by the finite-source effects. The black contours are derived from 
HST astrometry measurements assuming vTrei = 0.2 mas. The red contours show the joint 
probability distributions from both constraints. The lower two panels show the posterior 
probability distribution of the position angle 0^ of the relative lens-source proper motion 
for wide-binary and close-binary solutions, respectively. The histogram in red is derived 
from the red contours of joint probability for finite source and astrometry constraints in the 
upper panel. The blue histogram represents that of the microlens parallax. They mildly 
disagree at 2.5a. 
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Fig. 6. — Differences between OGLE V and HST F555W magnitudes for the matched stars 
are plotted against their V magnitudes measured by OGLE. To calculate the offset, we add 
a 0.017 mag "cosmic error" in quadrature to each point in order to reduce x^/dof to unity. 
The open circles represent the stars used to establish the final transformation, and the filled 
point shows an "outlier" . 
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Fig. 7. — Posterior probability distribution of lens mass M and relative lens-source paral- 
lax TTroi from MCMC simulations discussed in § 13.6. 1[ The constraints include those from 
parallax effects, finite-source effects and relative proper motion measurements from HST 
astrometry. The Ax^ = 1,4,9 contours are displayed in solid, dotted and dashed lines, 
respectively. Both wide-binary (magenta) and close-binary (blue) solutions are shown. The 
lines in black, red and green represent the predicted M and n^ev from the isochrones for 
different metal abundances: [M/H] = (black), —0.5 (red), —1.0 (green). The points on 
these lines correspond to the observed J-band magnitude / = 21.3 and various V — I values 
V — I = 1.8 (best estimate, filled dots), 2.0 (0.5 a, filled triangle), 2.1 (1.0 cr, filled squares), 
2.3 (1.5 cr, filled pentagons), and 2.6 (2.0 a, filled hexagons) 
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Fig. 8. — Posterior probability distribution of lens mass M and relative lens-source parallax 
TTrei from MCMC simulations assuming that the blended light comes from the lens star. The 
Ax^ = 1,4 contours are displayed in a solid hne for wide solutions, and in a dotted hne for 
close solutions. 




Fig. 9. — Probability contours of projected velocity (defined in Appendix E]) in the units 
of fc,± for both close-binary (upper panel) and wide-binary (lower panel) solutions. All the 
solutions that are outside the dotted circle are physically rejected as the velocities exceed the 
escape velocity of the system. The boundary in a solid line inside the dotted circle encloses 
the solutions for which circular orbits are allowed. 



- 46 - 




0.06 - 



^0.04 

CO 
O 



Ph 



0.02 - 







0.5 1 1.5 

log(a/AU) 

-| I I I I 11 I I I I I I I I I I I I I 



- Jn - I 



I I I I I I I I I I I I I I I I I I I 



20 40 60 80 100 

Temperature (K) 

1 I I I I I I I I I I I I I I I I I I I 




-1 -0.5 0.5 1 
cos(i) 



20 40 60 80 100 
Radial Velocity (m s~') 



Fig. 10. — Probability distributions of planetary parameters (semimajor axis a, equilibrium 
temperature, cosine of the inclination, and amplitude of radial velocity of the lens star) from 
MCMC realizations assuming circular orbital motion. Histograms in black and red represent 
the close-binary and wide-binary solutions, respectively. Dotted and dashed histograms 
represent the two degenerate solutions for each MCMC realization discussed in Appendix El 
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Fig. 11. — distributions for best-fit xallarap solutions at fixed binary-source orbital periods 
P. The solid and dotted lines represent xallarap fits with and without dynamical constraints 
described in § 13. 8[ The best-fit parallax solution is shown as a filled dot at period of 1 year. 
All of the fits shown in this figure assume no planetary orbital motion. 
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Fig. 12. — Results of xallarap fits by fixing binary orbital phase A and complement of 
inclination (3 at period P = 1 yr and Uq > 0. The plot is color-coded for solutions with Ax^ 
within 1 (black), 4 (red), 9 (green), 16 (blue), 25 (magenta), 49 (yellow) of the best fit. The 
Earth parameters are indicated by black circles. Because of a perfect symmetry {uq —uq 
and a — >■ —a), the upper black circle represents Earth parameter (A = 268°, (3 = —11°) for 
the case Uq < 0. Comparison of parallax with xallarap must be made with the better of the 
two, that is, the lower one. 



Table 1: Light Curve Parameter Estimations From Markov chain Monte Carlo Simulations. 



MCMC A 



Model 


to 






d 


9 


a 


P 






uo 


d/d 


Is 


/b 


Vs 






(HJD') 




(day) 




xlO^ 


(deg) 


XlO^ 










(mag) 


(mag) 


(mag) 


(mag) 


Wide+ 


3480.7024 


0.0282 


71.1 


1.306 


7.5 


273.63 


3.9 


-0.30 


-0.26 


-1.328 


-0.256 


19.51 


21.29 


20.85 


23.11 


1345.0 


+0.0058 
-0.0054 


+0.0008 
-0.0009 


+2.3 
-2.4 


+0.002 
-0.004 


±0.2 


+0.16 
-0.15 


+1.8 
-2.7 


+0.24 
-0.28 


±0.05 


+0.274 
-0.165 


+0.134 
-0.129 


+0.04 
-0.03 


+0.22 
-0.17 


±0.04 


+0.43 
-0.26 


Wide- 


3480.7028 


-0.0283 


70.6 


1.307 


7.5 


86.21 


3.9 


-0.34 


-0.26 


1.117 


-0.277 


19.51 


21.30 


20.85 


23.11 


1345.3 


+0.0054 
-0.0056 


+0.0010 
-0.0008 


±2.2 


+0.003 
-0.004 


±0.3 


+0.13 
-0.16 


+1.8 
-2.6 


+0.30 
-0.23 


±0.05 


+0.130 
-0.293 


+0.152 
-0.113 


±0.04 


±0.19 


±0.03 


+0.43 
-0.26 


Close+ 


3480.6789 


0.0239 


70.1 


0.763 


6.9 


274.27 


3.1 


-0.36 


-0.27 


0.301 


0.502 


19.52 


21.28 


20.85 


23.13 


1345.8 


+0.0055 
-0.0043 


+0.0009 
-0.0007 


+2.1 
-2.4 


+0.004 
-0.006 


±0.3 


+0.25 
-0.36 


+1.7 

-2.5 


+0.24 
-0.27 


±0.05 


+0.486 
-0.788 


+0.148 
-0.101 


±0.03 


+0.21 
-0.16 


+0.04 
-0.03 


+0.40 
-0.28 


Closc- 


3480.6799 


-0.0241 


69.2 


0.762 


6.9 


85.53 


2.7 


-0.33 


-0.26 


-0.405 


0.528 


19.52 


21.30 


20.85 


23.17 


1345.2 


+0.0042 
-0.0051 


+0.0009 
-0.0007 


+2.3 
-1.9 


+0.004 
-0.006 


±0.3 


+0.39 
-0.26 


±2.2 


+0.28 
-0.26 


±0.05 


+0.696 
-0.622 


+0.127 
-0.118 


±0.03 


±0.18 


±0.04 


+0.35 
-0.32 



MCMC B 



Model 


to 


Uq 


tE 


d 


Q 


a 


P 


tte.at 


tte.s 


u 


d/d 


Is 








x' 


(HJD') 




(day) 




xlO^ 


(deg) 


XlO^ 








(yr-^) 


(mag) 


(mag) 


(mag) 


(mag) 


Wide+ 


3480.7015 


0.0287 


69.3 


1.305 


7.7 


273.67 


6.1 


-0.02 


-0.22 


-1.242 


-0.283 


19.49 


21.40 


20.82 


23.91 


1353.4 


+0.0050 
-0.0059 


±0.0007 


+1.6 
-1.7 


+0.003 
-0.005 


±0.2 


+0.17 
-0.11 


±0.4 


±0.12 


±0.03 


+0.321 
-0.125 


±0.129 


±0.03 


±0.19 


±0.03 


+0.24 
-0.20 


Wide- 


3480.7012 


-0.0287 


69.2 


1.305 


7.7 


86.29 


6.0 


0.02 


-0.21 


1.193 


-0.293 


19.49 


21.40 


20.82 


23.97 


1353.3 


+0.0052 
-O.OOliO 


±0.0007 


+ 1.7 
-1.8 


+0.002 
-0.005 


±0.2 


+0.12 
-0.15 


+0.5 
-0.3 


+0.10 

-o.i:-! 


±0.03 


+0.127 
-0..S12 


+0,131 
-0.121 


±0.02 


±0.19 


±0.03 


+0,19 
-0.24 


Close± 


3480.6792 


0.0245 


68.3 


0.763 


7.0 


274.38 


6.0 


-0.01 


-0.22 


0.415 


0.569 


19.49 


21.35 


20.83 


23.88 


1355.5 


+0.0041 
-0.0051 


+0.0005 
-0.0006 


±1.6 


+0.003 
-0.006 


+0.3 
-0.2 


+0.23 
-0.39 


±0.4 


+0.12 
-0.15 


±0.02 


+0.503 
-0.744 


+0.112 
-0.130 


±0.03 


+0.20 
-0.16 


±0.02 


+0.24 

-0.18 


Closc- 


3480.6793 


-0.0245 


68.2 


0.762 


7.1 


85.63 


6.0 


0.04 


-0.22 


-0.179 


0.561 


19.50 


21.36 


20.83 


23.90 


1355.5 


+0.0042 
-0.0051 


±0.0006 


+1.8 
-1.5 


+0.004 
-0.006 


±0.3 


+0.46 
-0.24 


±0.4 


+0.09 
-0.16 


±0.03 


+0.703 
-0.722 


+0.126 
-0.112 


±0.02 


±0.18 


±0.02 


±0.21 
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Table 2: Derived Physical Parameters 



Model 


M 


TTrel 


Di 


/iN 


/iE 




Mp 






Me 


mas 


kpc 


mas yr^^ 


mas yr^"*^ 


mas 


^Jupiter 


AU 


Wide+ 


0.46 


0.19 


3.2 


-0.4 


-4.3 


0.84 


3.8 


3.6 


1353.4 


±0.04 


±0.04 


±0.4 


+2.7 
-3.1 


±0.3 


+0.06 
-0.04 


+0.3 
-0.4 


±0.2 


Wide- 


0.46 


0.19 


3.2 


0.3 


-4.3 


0.85 


3.8 


3.6 


1353.3 


±0.04 


+0.04 
-0.03 


±0.4 


+2.3 
-3.6 


+0.3 
-0.2 


±0.05 


+0.3 
-0.4 


±0.2 


Close+ 


0.46 


0.19 


3.1 


-2.6 


-4.4 


0.86 


3.4 


2.1 


1355.5 


±0.04 


+0.04 
-0.03 


±0.4 


+4.8 
-1.1 


±0.3 


±0.05 


+0.3 
-0.4 


±0.1 


Close— 


0.46 


0.20 


3.1 


-0.2 


-4.4 


0.87 


3.4 


2.1 


1355.5 


±0.04 


±0.04 


±0.3 


+3.6 
-3.4 


±0.3 


±0.04 


±0.3 


±0.1 
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Table 3: Jovian-mass Companions to M Dwarfs (M* < 0.55 Mq) 



Name 




Metallicity 


Dist. 


Mp 


P 


a 


Ref. 




(Mo) 




(pc) 


(Mj,p) 


(days) 


(AU) 




GJ 876c 


0.32 


-0.12 


4.660 


0.6- 


30.340 


0.13030 


1,2,3 




±0.03 


±0.12 


±0.004 


0.8 


±0.013 






GJ 876b 








1.9- 


60.940 


0.20783 












2.5 


±0.013 






GJ 849b 


0.49 


0.16 


8.8 


0.82/ sini 


1890 


2.35 


4 




±0.05 


±0.2 


±0.2 




±130 






GJ 317b 


0.24 


-0.23 


9.2 


1.2/ sini 


692.9 


0.95 


5 




±0.04 


±0.2 


±1.7 




±4 






GJ 832b 


0.45 


0.7 


~ 4.93 


0.64/ sini 


3416 


3.4 


6 




±0.05 


/-0.3 






±131 


±0.4 




OGLE-2006 


0.50 


? 


1490 


0.71 


1830 


2.3 


7 


-BLG-109Lb 


±0.05 


? 


±130 


±0.08 


±370 


±0.2 




OGLE-2006 








0.27 


5100 


4.6 




-BLG-lOQLc 








±0.03 


±730 


±0.5 




OGLE-2005 


0.46 


Subsolar?** 


3300 


3.8^ 




3.6^'" 


This 


-BLG-071Lb 


±0.04 




±300 


±0.4 




±0.2 


Paper 



"While the metaUicity of the OGLE-2005-BLG-071Lb host star is not directly constrained by our data, its 
kinematics indicate it is likely a member of the metal-poor thick disk. 

''We give the planet mass and projected separation for the wide solution, which is favored by Ax^ = 2.1. The 
second, close solution has Mp — 3.4 ± 0.3 Mjupitcr and r±^ = 2.1 ± 0.1 AU. 

'^We give the the projected separation between the host and planet at the time of event, which is the orbital 
parameter most directly constrained by our observations. However, assuming a circular orbit, we infer that 
the semi-major axis is likely only ~ 10 — 20% larger [a(wide) ~ 4.1AU, a(close) ~ 2.5AU]. 
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